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Abstract 

Discrete quantum mechanics is here defined to be a quantum theory of wave 
functions defined on integers Pi, Qi, while canonical quantum mechanics 
is assumed to be based on wave functions on the real numbers, R n . We 
study reversible mappings from the position operators and their quantum 
canonical operators pi of a canonical theory, onto the discrete, commuting 
operators Pi and Qi . In this paper we are particularly interested in harmonic 
oscillators. In the discrete system, these turn into deterministic models, which 
is our motivation for this study. We regard the procedure worked out here as 
a "canonical formalism" for discrete dynamics, and as a stepping stone to 
handling discrete deterministic systems in a quantum formalism. 
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1. Introduction 



In modern science, real numbers play such a fundamental role that it is difficult to imag- 
ine a world without real numbers. Nevertheless, one may suspect that real numbers are 
nothing but a human invention. By chance, humanity discovered over 2000 years ago that 
our world can be understood very accurately if we phraze its laws and its symmetries by 
manipulating real numbers, not only using addition and multiplication, but also subtrac- 
tion and division, and later of course also the extremely rich mathematical machinery 
beyond that, manipulations that do not work so well for integers alone, or even more 
limited quantities such as Boolean variables. 

Now imagine that, in contrast to these appearances, the real world, at its most fun- 
damental level, were not based on real numbers at all. We here consider systems where 
only the integers describe what happens at a deeper level. Can one understand why our 
world appears to be based on real numbers? 

The point we wish to make, and investigate, is that everything we customarily do with 
real numbers, can be done with integers also. A mapping exists that turns a set of two, 
mutually commuting, integer value operators, P and Q , into one real number valued 
operator. We can manipulate this real number anyway we like and subsequently map the 
result back onto the integers. However, the techniques to do this require the methods of 
quantum mechanics. Let us first give the general picture. 

Regardless what kind of theory we have, as long as its basic dymanical variables are 
defined to be sets of integers forming a space Z 2n , we can always assign an element of 
a basis of Hilbert space to each point in l? n . Having done that, the next thing we can 
do is transform to a different basis. For every single set of integers Q G Z , we can now 
consider the basis generated by the functions ip v (Q) — (Q\v) = e~ 27ra?< 2 , normalized by 
the orthogonality rule (771 1 772) = 5 (771 — 772) , if both 771 and 772 are chosen to lie within 
the unit interval (— |, |] . This is how any quantum theory on the set of integers can 
be mapped onto a quantum theory on the set of numbers 77 on the interval (— |, |] , or 
more precisely, on the unit circle. The mapping is unitary, so we can also have the inverse 
mapping. The hamiltonian on one side of the mapping is mapped onto a hamiltonian on 
the other side. 

Subsequently, one might consider deterministic theories. Since these would also either 
dictate how integers evolve, or how numbers evolve on a circle, these deterministic theories 
form a subclass of all quantum models on both sides of the mapping, but a deterministic 
model on one side in general will be mapped onto a pure quantum system on the other. 
Thus, our work may give further room for speculations about determinism in quantum 
mechanics, yet in this paper we are mainly interested in the mapping itself. 1 

How does the infinite line of real numbers emerge? Simply, we regard a real number 
as a composite of its integral part and a fractional part. The integral part, of course, is 
an integer, and the fractional part lives on the interval (— |, \ ] . This is the basis of our 
idea that wave functions defined on the real line can be transformed to functions on two 

1 A separate paper on deterministic versions of a quantum field theory is planned. 
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integers, P and Q , or two circles, rji and 7] 2 . If we know how we wish to transform these 
to a real coordinate q , we can calculate how its associated momentum p transforms, since 
p = —id/dq . All we will do now is work out the details, which, at first sight, appear to 
be straightforward. 

Intuitively, the physical situation may seem to be clear. If we put 2 [q,p] = i/27r , we 
have the "uncertainty relation" Sq ■ 5p ~ 1 ; we could take 8q ~ 1 and 5p ~ 1 , so that 
the integers Q and P could serve to enumerate a basis. If this is done with some care, 
the P and the Q can be made to commute. 

It so happened, however, that there are quite a bit of subtleties and pit falls, forcing 
us to go slow. To avoid making one of the many possible mistakes, we do our calculations 
slowly and carefully, which also exposes some beautiful underlying mathematics. 

We refer to this as a canonical formalism for discrete dynamics, complementary to 
the much more familiar canonical formalism of either classical mechanics or quantum 
mechanics. 



1.1. Notation 

While doing these calculations, the continuous need for factors 2n in our numerical 
expressions and normalization coefficients became irritating. When coordinates and mo- 
menta are treated symmetrically, factors y/2~K appear. A different, but somewhat unusual 
normalization can simplify things considerably. Instead of writing complex exponentials 
as e 2ma we will write 

e 2nia = e ia . £ = ^ _ 535 4Q .... ^ = \ if Z G Z . (1.1) 

Units of mass, time and energy will be normalized in such a way that 3 

2ir[x, p]=i, (x\e ipa \p) = (x + a\p) . (1.2) 

Throughout the paper, we will use 

capital Latin letters, N, P, Q, X, to indicate integers, (1.3) 

lower case Latin letters, p, q, x, to indicate real numbers, (1.4) 

and lower case Greek letters, a, r], £, A, for fractional numbers, (1-5) 

the latter being usually confined to the interval (— |, |] , but in expressions that are 
strictly periodic (with period 1 ), we may for simplicity replace this by the interval [0, 1) . 

States indicated as | P) and | Q) form a denumerable basis in Hilbert space. We in- 
terpret them as the momenta conjugated to the "position" operators operators | r] P ), \ tjq) 
on the periodic unit interval 4 : 

(v Q \Q)=e lQ ^, ( VP \P) = e iP ^ , (1.6) 



2 For the 2ir , see below under "notation". 

3 This means that, rather than h , it is Planck's original constant, h , that is normalized to one. 
4 Due to the unusual normalization of Eqs. (1.1), no factors containing 2tt are needed here. 
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(Qil Q2) = <W , rf Q ) = Stfg - rf Q ) , (1.7) 

and on the real line: (x\p) = e ipx , (xi\x 2 ) = S(xi — x 2 ) , (pi\P2) = ${pi — P2) ■ 

Furthermore, the operator t]q obeys: e lT) Q N \Q) = \Q+N) . On the unit interval (— |, |) , 
we can Fourier expand 

/\ j(—i) N 
i riidrn e- iNm = ^ K ^ e tNv , so that the t\q operator obeys 

The commutator between the operators i]q and Q is therefore not quite what one might 
have expected: 

(Qi\[vq, QWQ2) = (Q2 - Qi)(Qi\vq\Q2) = ± (<W - (-i)° 2 - Ql ) . (1.9) 

The extra term can be written as — ^IV'iX'^iI where the state l^i) , defined by (Q\ipi) — 
( — 1) Q , sits at the edge of the t]q interval. It arises because the t/q operator is not itself 
periodic, unless it is forced to jump back by one unit across the point i] = ±| , and this 
leads to a Dirac delta function in its derivative. We will refer to this state as an edge 
state. It will cause us quite some trouble, but, as we shall see, the edge states can be 
tamed. 



2. A first attempt 

Let us start doing what we intended to do from the very beginning: identify a position 
operator q as 

q = Q + Vp ; \q) = \Q,r)p). (2.1) 

This is totally legal, in principle; since P and Q are different integers, we work in the 
product Hilbert space, so that Q and r/p are commuting operators. All states \q) are 
represented exactly once, so that this identification is invertible. The question mark here 
indicates that we will replace this expression by a better choice later. 

Computation of the associated \p) states is straightforward: 

1 

/co poo ? 

dq\q)(q\p) = / dqe^\q) = ^ 1 ^ p(Q+Vp) \Q,Vp) ■ (2.2) 
-00 J — 00 q J — ~2 

Let us write p — K + k , where, in line with the notation that we will always use, K is 
an integer and \k\ < \. Now we have the general integral expression 

/\ 2 sin 7tk ( — 1 ) N 
drje i{N+K)r > = — , N integer, « fractional, (2.3) 

'2 
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(which, In the limit k — > , turns into 5^ ), and from this we find 

while from (2.1), we have: (Qi,Pi\Q + Vp) = <We~ ir?pPl • (2-5 

For future use, we also need the corresponding expressions in r\Q,r\ P space. One easily 
derives 

(tj^tjUK + k) I StfQ + K)^**"). (2.6 

(vhMQ + O = s( Vp -0^ Q . (2.7 

A drawback may seem to be that our treatment of the operators p and q is asym- 
metric, but this is easy to remedy; we can replace Eqs. (2.6) and (2.7) by 

(v^VpIK + k) = S{ V ^ + K)^ K+m \ (2.8 

(vh>1p\Q + t) = ^ P -i)^ {Q+n \ (2.9 
and verify that the first parts of Eq. (2.2) still hold if the parameters a and a obey 

a + a=l. (2.10 

We obtain 

(0,P,\K + k) - 2_sin^J-l)^ 

(o P\o + e\ - 2sin (^) (-i) Q - Ql (212 

and symmetry is obtained when a — a — \ . 

One might have hoped for a formalism where K coincides with P 1 and Q with Qi , 
but we see that this cannot be made generally true. The kernels, which vanish no faster 
than 1/\K — P\\ or 1/\Q — Q\ \ , give problems when we wish to construct the operators 
p and q themselves: non-canonical terms appear in their commutator. Again, all this is 
due to edge states. One way to avoid disaster is by only using soft operators such as e zap 
and e %l3q , to be constructed by using 

i 

(Qi,Pi\e iap \Q 2 ,P 2 ) = Y, f\_MQiiPi\K + K)^ K+K \K+K\Q 2 ,P 2 ) , (2.13) 

K "2 

and so on, since these sums and integrals converge well. Indeed, one then obtains the 
desired commutation rule 

^iap e if3x _ e ia(3 ^iap ^ 14) 

This however does not allow us to compute hamiltonians and other interesting operators 
directly in (P, Q) space, which is why a superior procedure was searched for. It turns 
out that all of the above definitions for the p and q states can be replaced by better 
expressions, which is why we adorned them with question marks. 
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3. (Almost) removing the edge states 

If the p and q states of the previous section are used to compute the commutator [q, p] , 
one finds that, at the edges of the intervals where r/Q and r/p are defined, the edge states 
give non-canonical contributions. Although there is nothing wrong with defining our 
mappings this way, it turns out that the contributions of the edge states are cumbersome. 
It is better to reduce their effects to a minimum. To achieve this, we have to recover full 
periodicity in i]q and r]p . For technical reasons, we rename these as 771 and i] 2 . 

Consider the torus defined by the intervals — \ < 771 < \ and — | < 772 < §. In 
order to avoid all contributions from edge states, we need to define the operators q and 
p unambiguously on this torus, with full periodicity. From the expressions (2.8)-(2.12), 
we are led to try 5 : 

q = -^d 1 + ar] 2 , p=-^d 2 -ar) 1 \ d 1 ^d/dr] 1 , d 2 = d/dr} 2 , (3.1) 

with a + a = 1 . Indeed, since 

[d^rjj] = 8 i:i , [771, 772] = , and d 2 ] = , (3.2) 

we find that Eq. (3.1) ensures that [q, p] — ^- . However, if we also desire that q = 
Q + 772 and p — K — r\\ , where Q and K are always integral then the wave functions 
ipq P (r]i,r) 2 ) allowed on the torus must always obey the following continuity properties 
across the borders: 

<M + i V2) = V>qp(-i V2)e^ 2 ; Vqp(?7l, +|) = ^qpfal, -^e^ . (3.3) 

The q and p eigenstates are then: 

If q = X + i : (77 1 , 77 2 |g) 
If p = K + k : (77i, 77 2 |p) 

Now, observe something that is of crucial importance: the only way to introduce on 
the torus the operators Q = —idi and P = —id 2 , where both Q and P are restricted 
to be integers, is to have regular periodicity. We need a Hilbert space of wave functions 
^PQ (?7i,?72) that obey: 

#q(+|, m) = ^pq(-i V2) ; 4>pq(Vi, +|) = 4>pq(Vi, ~\) ■ (3-6) 

To relate the functions ip qp to -^pq , we can either introduce a factor e ta ' ni ' n2 , so that 
the first part of Eq. (3.3) agrees with the first part of Eq. (3.6), or a factor e~' l(Tr]xrn to 
make the seconds parts agree, but to make them all agree is harder. The interpolating 
function on the torus must be unimodular if we want the transformation to be unitary, 
while its phase must make a full rotation over 2n while circulating around the edge 
if the unit square — in short: the interpolating factor must contain one unit of phase 
flux. This gives the interpolating function a singularity, which is easiest to situate on the 



= $(V2 - i)e iX ^ +iar '^ , and (3.4) 
= 5U + K ) e iKr *- iar >w . (3.5) 
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Figure 1: An artist's impression of the contour lines of the function 0(?7,£) , defined in 
Eq. (3.14). Its values run from —j (top left) to +| (top right), corresponding to the 
angles -90° to 90°. 



corner: (rji — ±~, 772 = ±§) ■ The required properties of this phase function e 1 ^ 1 ^ are 
illustrated in Fig. 1. 

With such an interpolating phase function, we can map the real line on the double 
set of integers (Q, P) . The fact that there is a singularity means that there is one edge 
state left: 

fal.^edge) = &(m - §)<% 2 - \) , (Q, P |^ edge ) = (-lf- Q • (3.7) 

This has to be contrasted with the situation we had before: there was an edge state 
{Q l^i) = ( — 1) Q f° r au values of P and an edge state {P — (~ 1) P f° r ah values of 
Q . We now have just a single edge state left. This improves convergence sufficiently to 
allow us to move on and study the q and p states in the (Q, P) Hilbert space without 
further trouble. The single edge state cannot be avoided since it is associated with a 
conserved flux on the torus. 

We do have to construct explicit expressions for the interpolating function. The desired 
function is constructed in two steps. First, find a smooth complex function / 772) 
obeying both boundary condtions (3.3). A suitable choice could be: 

/(771, Tfe) = e~ iWmm cosTTT/2 + e ia ™ cosvrr/! , (3.8) 

and we can confirm here that this function can serve our purpose quite well. It can 
be used to connect functions with the matching conditions (3.3) with functions obeying 

5 Later we will observe that no generality is lost if we return to the case £7 = 1, £7 = 0, but keeping 
the two parameters cr, a makes it easier to observe the p f-> q symmetry 
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(3.6). However, it is not the optimal choice, since the border crossings are not completely 
analytic. We found a better solution in terms of the elliptic theta functions [1]. We choose: 



i 



f( m ,r] 2 ) = e"™ 2 e ~* N +N{ivi+m) ■ (3.9) 

N=-oo 

Within the intervals |?7i| < \ and 1^1 < \ , the functions (3.8) and (3.9) are so similar 
that their effects are very nearly the same, but they do not coincide. 

The interpolating function is then 

u(m,V2) = f}? uri2 l, ; ^pq{vuV2) = ^(^1,^2)^(^1,^2) • (3. 10) 

Note that, if 771 = \ +a and r] 2 = \ +b, with a and b infinitesimal, then e ™vim/ e -i°mm = 
i, while cos7r^! = — na and cos7T77 2 = —nb, so the function / in Eq. (3.8) makes a full 
phase rotation if we follow a small curve around the point (|, |) . The function in Eq. (3.9) 
does the same (See Appendix A). Both functions have a zero in the corners (±^, ±|) • 

The operators {i]Q,i]p) are now defined to take the same values as (771,772) , but we 
define their joint fundamental eigen states to be rotated by a factor f/ (771 , 772) • Thus, we 
rewrite (3.10) as follows: 

\Vq,Vp) = U*^, V2) \VuV2) ■ (3.11) 
Thus, Eqs. (3.4) and (3.5) for the p and q eigen states are now written as 

(VQ,Vp\q) = U(r ]Q ,r ]P )S(r ]P -0e iXl,Q+iWriQr ' p , g = A + £; and (3.12) 
(VQ,Vp\p) = U(r ]Q ,r ]P )5(r ]Q + Ky K ^- i ^ , p = K + k . (3.13) 

Note that a (a) dependence in Eqs. (3.12), (3.13), (3.8) and (3.9) cancels out. We had 
kept a and a just to demonstrate the logical coherence of these equations. Clearly, the 
p -H- q symmetry is fully maintained in our present procedure. 

Writing 6 : 

j2 e -^ + N^ H) = rM € i<P(v,0 > (3 . 14) 

N 

with r and real, we need the phase function 0(77, £). Its most important properties 
are that it is differentiable, and, when using (3.9), we also have 

0(r7, £ + 1) = 0(77, + 77; 0(77 + 1, = 0(77, ; (3.15) 

0(77, O = -0H7,O = -0(77,-0 • (3.16) 



6 We apologise for the use of brackets ( ) that can mean two different things; most often they are just 
meant to group terms together when multiplied, but in expressions such as (f)(r], £) they indicate that <j> 
is a function of 77 and £ . The comma should make this unambiguous. 
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One can also prove the important relation: 

<i>(v,0 + <f>(Z,v) = Zv, (3.17) 

and in Eq. (3.11) we use: 

U(t]i,t]2) = e iffrnm ~ i ^ Vl ' m ^ = e ~ ^0 " ^ ' 1 '? 2 + ^< ^ ^, ( , ? 2 > , ? 1 ) _ (3.18) 

Eqs. (3.12) and (3.13) have now become 
(VQ, VM = e iX ^ +i ^^ 5(rj P - ; (77o, Vp\p) = e iK ^^ S(rj Q + k) . (3.19) 

Since, in these equations, the states \q) and \p) are now fully periodic in r\ P and i]q , it 
finally became legitimate to expand in the basis \Q,P) = e' lPr,p+ ' lQ ' nQ \r]Q,rip} , to obtain 

i i 
(Q, P\q) = e- ip t j\ drie i<Kv,0+i(x-Q)v . (Q ? P \ p ) = e <Q« d?7 e *<K^)+^-^ _ (3>20) 

"2 "2 

This expression for (0,0\q) can be inverted to give 

e iHv,0 = ^(o, 0|X + e"**" . (3.21) 

To illustrate what was gained by this new mapping between \q) and \p) states on the 
one hand and the \Q,P) states on the other, the wave functions (g|0,0) and (p\0,0) in 
various schemes are illustrated in Fig. 2. The Figure suggests that the function (0,0|g) , 
for large values of \q\ , tends to 

^ (-l) X 5(\q\ -X- I) , x 

£ WTiy (3 ' 22) 

The matrix elements (Q,P\q) and (Q,P\p) can all be easily expressed in terms of the 
one basic function ip(q) = (0, 0\q) . Further properties of this function are expanded upon 
in Appendix A. 

The contour lines of the function 0(t?, £) are sketched in Fig. I. Observe the singularity 
due to the edge state in the corners. In (r)Q,r)p) space, the q and p-operators now read: 

9 = ^4K4' A(,p ' ,a) ) = S^ + " p -(^ ( "°'" p) ) ; (X23) 
p = ^Wp'iwp^-^) = ^9^"" e + (^ ( '' p ''' o) ) ' (3 - 24) 

where we used Eq. (3.17). Note that, according to Eqs. (3.16), these expressions are 
exactly periodic in both tjq and rjp , with the only singularities being in the corners 
of the (t)q, rjp) -quadrant. Apart from possible effects due to this singularity, we have, 
everywhere in (r)Q,r)p) space, 

[q,p}=t/2n . (3.25) 

What was gained by adding the <fi field is, that there are no contributions from 5 functions 
on the boundaries of our quadrant. Nevetheless, the effects of the one remaining edge 
state still show up in the commutator (3.25) when expressed in the Q, P basis. The 
relevant matrix elements are computed in Appendix B. 
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a) 




Figure 2: \Q,P) states in various cases, a) Our starting point, he state 1 , 0) in q 
space, Eq. (2.5) . b) Its Fourier transform, in p space, Eq. (2.4), converges only slowly 
for large p. c) The 1 , 0) state according to Eq. (2.12). It is equal to its own Fourier 
transform (its form in p space), d) The 1 , 0) state according to Eq. (3.20). It is also 
equal to its Fourier transform in p space, and converges better. The other \Q,P) states 
are obtained from these by translation and/or multiplication with e lPq . e) The small 
peaks in the curve d quickly reach this asymptotic form, if multiplied both horizontally 
and vertically with a factor (X + .5) . The surface areas of these peaks rapidly approach 
the values ±1/(2tt(X + |) 2 ) . 

4. The harmonic oscillator 

We are interested in the fate of the quantum harmonic oscillator in p-q space when 
we map it to the discrete P, Q variables. To find out, we first have to go the f]Q-,r\p 
quadrant. In our units, the hamiltonian is 

H = n(p 2 + q 2 ) = a)a + \ ; |V = -2iriHif> , i/i(t) = e~ iHt ^j(0) . (4.1) 

Creation and annihilation operators are 

cr = Vtt(p + iq) , a = \fn(p — iq) ; [a, cr] = 1 . (4.2) 

The eigen states \ip n ) , n = 0, 1, • • • , are given by 

a|^o)=0, 1^) = ^"!^). (4.3) 



9 



In terms of t/q and 7] P , for simplicity again written as 771,772 , the annihilation operator 



is 



a = V^(-^(di+id 2 )-i V2 + ±(d 1 + id 2 ) ( f ) (rj 1: 772)) . (4.4) 

Solving the differential equation a^o) = appears to be easy. Introduce z± = 
771 ± 7772 , to find 

cL(lV>o)e^) = id-<f>( m , m ) |^o) ; fai,^) = /(z + )e^ 2) -^ 2 . (4.5) 

The function f(z+) = /(771 + 2772) is now determined by the boundary conditions (3.16). 
The term |7/| plays no role there, because it matches periodic boundary conditions, but 
the behavior of <fi is more difficult to accommodate for. Fortunately, of course, we know 
the solution in q space: 

<<#o> = 2 1/4 e"^ 2 , (4.6) 

So, using Eq. (3.19), 

(77!, 772^0) = 2 1 / 4 Ld^-^e-^^^^ = (4.7) 

X=-oo 2 



which is indeed of the form (4.5). 

The expression (4.8) could have been arrived at directly by the following chain of 
arguments: the periodicity requirements (3.16) imply strict periodicity in r/i ; therefore, 
the wave function can be expanded in waves periodic in r/i : }\z + ) = Xlxez a x (V2) ^ Xr>1 ■ 
This must be a function of 771 + 7772 ; therefore, 

f( z+ ) = J2^x^ X{vi+lV2) ■ (4-9) 
x 

According to the first equation (3.16), we must require 

f(Vi + ¥) = f(Vi ~ 5 therefore, (4.10) 

1 „ * _ , 1 . . _ .,„, „. ,1 . 



J2 a x£ Xm ~i x = J2 a x etXm+ * x ~ 1V1 = J2 a x ei{x ~ 1)m+ * x ( 4 - n ) 

XX X 

= + or (4.12) 



x 



a x e~ X -i ; a x = Ce~2 x \ (4.13) 



l x + 1 
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To construct the Q, P matrix elements of the hamiltonian (4.1), it seems to be best 
to use the wave functions of Fig. 2 d , 

(q\Q, P) = e iP « (q - Q\0, 0) and (p\Q, P) = e"** (p - P\0, 0) , (4.14) 

as templates — note, that the functions (q\0,0) and (p|0,0) are given by the same 
mathematical expression. Then the matrix elements of H = n(q 2 + p 2 ) are found by 
calculating 

/oo 
dqe iq ^- p ^(0,0\q - Q 1 )q 2 (q - Q 2 \0,0) ; (4.15) 
■oo 

/oo 
dpe^ 1 -^)(0,0b-P 1 )p 2 (p-P 2 |0,0) , (4.16) 
-oo 

where the curves of Fig. 2 must be inserted. Considering the asymptotic form (3.22), we 
notice that both integrals here diverge logarithmically. This requires a logarithmically 
infinite subtraction in these matrix elements, proportional to (— l^i-Qa+Pi-ft ^ w hich 
is exactly the contribution of the one remaining edge state, Eq. (3.7). It is a harmless 
additive constant in H . 

The hamiltonian (4.1) describes a harmonic oscillator with period T = 1 , but in terms 
of the discrete states \P, Q) , the evolution is deterministic over multiples of one quarter 
of this period, as will be demonstrated now. As is well-known for a harmonic oscilator, it 
sends position states into momentum states and vice versa 7 . If at t — we have a state 

<fl#M>) = /(?) , (4-17) 
then at one quarter of a period one finds 

(*#t=i/4) = f(-p) ■ (4.18) 

and so on. So now, because of the symmetry in the states \p) and \q) in the expressions 
(4.14), we find that if the initial state is 

(Q,P\A=o) = 5 QA 5 PB , (4.19) 

then after one quarter of a period, 

(Q,P|^=i/4> = Wp,-a, ( 4 - 2 °) 

and so on. This is a kind of determinism that we plan to study further: Q turns into P 
and P turns into — Q . What is new in our formalism is that we identified a hamiltonian 
that does this job in Q, P space, while it has a natural and non-trivial ground state, 
Eq. (4.8). 

7 In fact, the harmonic oscillator is the ideal instrument to produce fractional Fourier trans forms [2], 
by considering how a wave function transforms at arbitrary, fractional time t . 
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5. The edge state 



The one edge state of our system is ^> e d g e , described by 

(VQ,Vp\^e dg e) = 6(r) Q - \) 6(r) P - \) or (Q, Pledge) = (-1) P+Q . (5.1) 

Here, we did not attempt to normalize it; the norm squared of this state would be LqL p , if 
we would have had a box of length Lq in the Q direction and Lp in the P direction. It is 
this state that causes a logarithmic divergence in the hamiltonian of a harmanic oscillator, 
but more importantly, it generates a non-canonical term in the commutation rule for the 
real number operators q and p that we constructed. In Eq. (B.13) in Appendix B, we 
see that the non-canonical contribution to this commutator is proportional to 

|^edge)(^edge| , (5.2) 



Thus, to achieve a useful mapping with a standard canonical system, somehow, this 
state must be filtered out. It would be tempting to argue that we should simply replace 
all states by 

IVWs) = ( 1 - Z^LF l^edgeXV'edgel )l/>) , (5-3) 

but this is not the entire story. When acting with operators such as q and p on such a 
"physical" state, a non-physical state might emerge. A non-physical state is a state 
that obeys 

^(-l) Q+p (Q,P|^a(V0^O. (5.4) 

P,Q 

There is a much better way to look at the edge state. To see how to remove an 
unphysical state correctly, we now formulate an important theorem: 

Let be a state with the following properties: 

1. It has a compact support in (Q, P) space, or equivalently, for only a finite 
number of values Q and P , its components (Q, P\ip) ^ , and 

2. Its "edge state coefficient" a{jp) , as defined in Eq. (5.4), is not equal to 
zero, 

then the operator p 2 acting on this state has a positive, infinitely large ex- 
pectation value: {ip\p 2 \ip) — > oo . 

To prove this theorem, just consider the matrix elements (B.10), (B.ll) of p, derived in 
Appendix B, and use them to derive the asymptotic form of the matrix elements 

(Q,p\ P m = J2 (Q,p\p\Q2,P2){Q2,P2m , (5.5) 

Q 2 ,P 2 



12 



for large values of Q and P . Because all values of Q2 and P2 are bounded, we find that 
only the "vector potential" operator ap contributes for large Q and P , and these Q,P 
matrix elements approach the values 

t^w • (5 - 6) 

Therefore, the expectation value of p 2 receives the asymptotic contribution 



{ h( p2 + Q 



2 



(5.7) 



and the sum of these positive terms diverges logarithmically. Q.E.D. 

Now, we use this result as follows. If we decide to map our discrete system on a 
continuous model where the hamiltonian has the form 

H^^ P 2 + V(q), (5.8) 

where V(q) is bounded from below, then any state ip with a(tp) 7^ 0, has an infinite 
expectation value of this hamiltonian. If we would only be interested in all those states 
for which the hamiltonian is finite, these will automatically obey a(ip) = . So, for these 
"physical" states, we can also ignore the non-canonical term in the [q,p] commutator 
(the -1 in Eq. (B.13)). 

We now see why also our harmonic oscillator has a divergent term in its hamiltonian, 
Eqs. (4.15), (4.16), proportional to a(jp) . It is not difficult to generalize our result for 
hamiltonians that contain a magnetic term: H — ^(p — A(q)) 2 + V(q) . 

Clearly, all physical states (those states that have finite energy in the canonical (q,p) 
system) have the property that, in (tjq, rjp) space, they vanish at the points in the corners: 
t)q — ±|, rjp = ±| . This is a non-trivial property; it is at these corners where the vector 
potential field (a,Q,ap) has a singularity. 

Notice furthermore, that our argument would not be so straightforward if we had edge 
states for all values of either P or Q (that is, a singularity at the entire boundary of 
{VQiVp) space, instead of only at the corners), which is what our first attempts were 
giving us (see Eq. (1.9)). It is much better to have just one edge state. 



6. Discussion 

The study carried out in this paper could be viewed as a simple exercise in the application 
of matrices to transform from one basis of Hilbert space to another. However, it was done 
with a very special purpose, which is to demonstrate that quantum theories based on real 
numbers, such as the popsitions and momenta of elementary particles, or equivalently, the 
real values of bosonic field variables, can be mapped onto theories where the fundamental 
degrees of freedom are integers. Our aim is to apply this to cellular automaton theories 
of Nature, possibly of importance at or near the Planck scale. These theories would first 
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have to be cast in the form of one of the more familiar quantum field theories by means 
of our transformation matrices, after which the renormalization group can be applied to 
determine their behavior at very large distance and low energy scales. 

Of course, one would primarily think of quantum theories acting on integer-valued 
oprators, but the author suspects that one may even go one (important) step further. The 
deterministic evolution of a classical cellular automaton may be cast in a mathematical 
framework where Hilbert space is introduced artificially [3], just to describe the evolution 
process as time proceeds. Only at integral time steps, this evolution would take the form 
of pure permutations in our discrete basis elements; due to the technique of introducing 
Hilbert space however, one can turn the time variable into a continuous variable of the kind 
that we think we experience in the physical world. Allowing ourselves the complete set of 
Hilbert space transformations to any basis we like, makes these theories mathematically 
hardly distinguishable from conventional quantum mechanics. Indeed, after applying 
these techniques, distinguishing deterministic theories from quantum mechanics may be 
so difficult at larger scales that this could be an explanation as to why we think our world 
is quantum mechanical, while it may actually be deterministic at the Planck scale. 

We are aware of the fact that such arguments seem to be completely irreconcilable 
with the numerous examples of experimentally realizable situations[4] where Bell's in- 
equalities^] are violated. To analyze exhaustively how this apparent disagreement can 
be resolved is beyond the aim of this paper (see for instance Ref. [6]). We would just 
emphasize the following facts: 

- Our quantum states are real quantum states, including Born's identification of am- 
plitudes squared with probabilities. There is nothing against the use of states of 
this kind that are as entangled as in real or imagined Aspect-like experiments, even 
if the underlying theory happens to be a deterministic one. 

- When considering Bell's inequalities in a deterministic system, in the real world, this 
system cannot be in just any quantum state; it is in exactly one of the fundamental 
cellular automaton states. In the real world no superpositions of such states can 
occur. 

- As stated however, a complete analysis is not possible without more detailed models. 
These we do not show here, but more will come in our next paper. 

In Section 4, a brief illustration was given of a simple deterministic automaton. This 
particular example turns back to its original position after just four steps. The importance 
of this model is not the demonstration that it exists, since such models are quite trivial; 
the hamiltonian H merely needs to obey the condition that 

e - 4iH = i ; h = \N + integer , (6.1) 

where the value of the integer is formally immaterial. 

The importance of this model is that we gave a very special representation of the 
hamiltonian, where the integer takes values such that we have a real harmonic oscillator. 
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In particular, there is a natural ground state. This is important because it allows us to do 
thermodynamics with this oscillator. In our previous work, the choice of a good hamilto- 
nian, with a non-trivial and physically meaningful ground state, has been problematic. 

Now what we have is a mapping of models defined in a discrete world, onto models 
that are quantum mechanical in the conventional sense, with pi and g« operators obeying 
[QiiPj] = ifiij ■ It sould be emphasized that there are different ways to perform this 
mapping, some of which may seem to be easier mathematically and conceptually, but we 
think that the procedure we ended up with is unique and superior. It is totally symmetric 
in p and q , it has only one edge state, and, in spite of the rather awkward mathematical 
expression for the template states (0,0\q) (Appendix A), we ended up in remarkably 
simple expressions for the q and p matrix elements in (Q,P) space (see Appendix B), 
so that the whole scheme ends up in being quite transparent. 

Now our real aim is to apply this mapping in the case that the model in (Q, P) space 
is a deterministic one. This means that, at integral time intervals t (if the fundamental 
time unit is normalized to 1), the evolution U{t) is a pure permutation (possibly a quite 
complicated one) on the (Q, P) lattice. Although the observables are only well-defined 
and "ontological" at integral time intervals, it is relatively easy to devise a hamiltonian 
H such that 



This hamiltonian then, allows us to extrapolate the time variable to be continuous, but 
we think it is reasonable to postulate that, if the time quantum is sufficiently rapid, a 
"physical observer" cannot distinguish between integer time states and non-integer time 
states. 

It is an essential feature that, if we would add arbitrary integers to any of the energy 
eigenvalues, the original discrete model will not be affected at all, whereas the canonical 
q,p model may seem to become totally different. This means that, although our mapping 
{Q,P) -H- (q,p) is unique and inversible, there are many canonical hamiltonians that all 
correspond to the same discrete system. This freedom may actually be used to seek for 
a hamiltonian that is as close as possible to one of the standard quantum mechanical 
systems, such as a quantum field theory. 

Work is under way on a sequel of this paper, where we construct such a physically much 
more interesting model: a quantum field theory in one space- and one time dimension, 
having a non-trivial hamiltonian with a ground state; we suspect that the significance of 
the present procedures will then become more evident than it may be now. 

A. Properties of the function (0,0\q) . 



To go from the real numbers q and p to the integers Q and P , we use the real function 
defined in Eqs. (3.20) and (3.21), 



U(t) = e 



-iHt 




#z)=^(-<z) = <0,0|g) = <g|0,0), 



(A.I) 
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which is equal to its own Fourier transform, 

${p) = J dq(q\0,0)e- ipq . (A.2) 

The other matrix elements are simply given by 

(Q,P\q) = ii(q-Q)e- iP « . (A.3) 

ip(q) is sketched in Fig. 2, and it plays a central role in our mappings. Because of the 
periodicity properties (3.15), the definition of the function ip can be written as 

1 

Tj){x) = drie l ^ x) , (A.4) 

and the function (j)(r),x) is given by Eq. (3.14), or 

00 i 

r(r j ,x)e i ^ r '' x) = e-2 K2+K ( x+ir » ; r, real. (A.5) 

K=-oo 

This sum is a special case of the elliptic function $3 , and it can also be written as a 
product: 

00 00 

r{i 1) x)e l<t> ^ x) = ]J(l-e- K ) f[ (1 + e x+vr, ~ K ~2 )(1 + e ~ x ~ 2) , (A.6) 
AT=1 x=o 

with r and real. Here, the first product term is of lesser importance since it only 
multiplies r(r),x) with a constant, while not contributing to (/>(?], x) . Note, that <f)(i],x) 
has a vortex singularity when e 1 ^^ has a zero, and these zeros can easily be read off from 
Eq. (A.6); they are located at (rj, x) = (K 1 + \, K 2 + \) ■ We see that in Eq. (A.4), the 
absolute value r(r),x) of the sum in (A.5), or the product in (A.6), has been divided out, 
and this makes the evaluation of the integral over rj hard, although it is well bounded. 

In Eq. (A.5), the sum is dominated by the K value closest to x . In Fig. 2, the small 
peaks at large x (Fig. 2 e ) arise when the dominant K value in the sum switches from 
one integer to the next. 

Unitarity property : 

From the fact that e l ^ v ' x ^ in Eq. (A.4) is the Fourier transform of ip( x ) f° r integral p, 
and that it has absolute value one, we derive that 

00 

^2 ip{x + K)ip(x + K + M) = 5 M o , K,Lez (A.7) 

K=-oo 

(use was made of Eq. (3.21)). 
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B. Matrix elements of q and p operators 



The matrix elements (Q±, P±\ q \Q 2 , P 2 ) can be calculated explicitly. Let us first compute 
the operator q in tjq, r)p -space. We write the expression (3.23) as follows: 

-* d , / \ , \ d(f)(r}p,r} Q ) 



2tt dr] Q dr] Q 

where the function oq is regarded as the Q -component of a vector potential field a . For 
the phase 4>{f)p,r]Q) , we can now best use the product formula (A. 6), which gives: 

00 

(Iq(vq,Vp) = ^2^(r] Q ,r]p) , (B.2) 

K=0 

a Q<M>Vp) = (arg(l + 6^+^-1) + arg(l + e-^-^-^l)) . (B.3) 



Evaluation gives: 
aq(VQ,Vp) 



|sin(27r?7p) , | sm(27irj P ) 



cos(27r?7p) + cosh(27r(r]Q — X — |)) cos(27t?7p) + cosh(27r(r]Q + K + |)) 

which can now be rewritten in a more compact way by rewriting Eq. (B.2) as a sum for 
K values running from —00 to 00 instead of to 00 . 

By writing 

(P 1 \r ]P }(r ]P \P 2 } = e tP , P = P 2 - P 1 , (B.5) 

we now proceed to write the matrix elements of the operator ag in the (tjq, P) frame: 

00 

(P 1 |a Q (7 7Q )|P 2 ) = a Q^ P ' K ) ' P = P2~Pi, (B.6) 

K=-oo 

finding 



( VQ , P, K) = \sgn(P)(-l) p -He~\ p ^ +K+ ^\ , (B.7) 



where sgn(P) is defined to be ±1 if P < and if P = . The absolute value taken in 
the exponent indeed means that we always have a negative exponent there; it originated 
when the contour integral forced us to choose a pole inside the unit circle. 

Next, we find the (Q, P) matrix elements by integrating this with a factor e tQ , with 
Q — Q2 — Qi , to obtain the remarkably simple expression 

(-i) p +Q+ 1 iP 

(Qi,Pi|oq|Q2,P2> = 1 p 2 + Q2 ■ (B.8) 
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In Eq. (B.l) this gives for the q operator: 

q = Q + a Q ; (Q u Pi\q\Q 2 , P 2 ) = Qi<W S PlP2 + (Q u P 1 \a Q \Q 2 , P 2 ) . (B.9) 

For the p operator, one obtains analogously, writing P = P 2 — Pi , 

p = P + a P , (B.10) 
(Qi,Pi|ap|Q 2 ,P 2 > = p 2 J + g2 ■ (B.11) 

It is important to check the commutation rule for g and p . Doing the matrix multi- 
plications for the matrices (B.9) and (B.10), one finds that 

[Q,P} = 0, [a Q ,a P ] = 0, (B.12) 
(Q 1 ,P 1 \[q,p]\Q 2 ,P 2 ) = (Q 1 ,P 1 \[Q,a P ] + [a Q ,P]\Q 2 ,P 2 ) = 

i(-l) Q i- Q > +p i- p *(5 QlQ2 5 Pl P 2 ~ !)• (B.13) 

Again, we see that the desired commutation rule, [q, p] — % , is obeyed only after we project 
out the edge state ?/>edge by demanding that all our states must obey (V'edgel i/)) = , see 
Eq. (3.7), see Section 5. 
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